Near-threshold production of the multi-strange E hyperon 
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The yield for the multi-strange H~ hyperon has been measured in 6 AGeV Au+Au collisions 
via reconstruction of its decay products vr" and A, the latter also being reconstructed from its 
daughter tracks of tt~ and p. The measurement is rather close to the threshold for S~ production 
and therefore provides an important test of model predictions. The measured yield for S~ and A 
are compared for several centralities. In central collisions the H~ yield is found to be in excellent 
agreement with statistical and transport model predictions, suggesting that multi-strange hadron 
production approaches chemical equilibrium in high baryon density nuclear matter. 
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Many years ago, it was proposed that strange particle 
yields could serve as an indicator for the formation of the 
quark-gluon plasma (QGP) [1]. For a partonic system 
with partially restored chiral symmetry, strangeness is 
expected to be readily produced due to the lower energy 
threshold for ss quark pair production. The production 
of strangeness may be further increased for significant 
baryon densities where the chemical potential associated 
with the production of light quarks is raised. For these 
reasons it is expected that the QGP phase should re- 
flect a high strange (anti)quark density [2]. Due to their 
strange quark content, multi-strange baryons [and anti- 
baryons] are expected to be a more sensitive probe of 
this phase than hadrons which contain only one strange 
valence quark. 

Interest in strangeness and especially in the yield of 
multi-strange hadrons has grown enormously [3]. The 
observation of "strangeness enhancement" in the produc- 
tion of S~ (dss) and fl~ {sss) in heavy ion collisions 



compared to pp or pA collisions has contributed to the 
excitement [4]. 

A comparison of experimental and calculated yields 
indicates two distinct features: (i) Statistical models, 
with input of temperature and chemical potential, are 
in fair or good agreement with the yields reported from 
the Alternating Gradient Synchrotron (AGS) [5], the 
Super Proton Synchrotron (SPS) [6] and from recent 
experiments performed at the Relativistic Heavy Ion 
Collider (RHIC) [7,8]. Using a canonical description 
this approach is successful even for beam energies of 
~ 1 AGeV [9]. They can even explain the observed 
enhancement between pA and Pb+Pb collisions [10]. 
(ii) Transport-model calculations based on hadronic in- 
teractions are currently unable to account for the mea- 
sured yields of multi-strange particles at SPS energies 
(158 AGeV) [11]. The success of the statistical model is 
surprising and raises the question as to whether or not 
this success holds true for near-threshold multi-strange 
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particle production. 

This letter reports on the production of the singly- 
strange A hyperon and the doubly-strange S~ hyperon in 
6 AGcV Au+Au collisions. Preliminary reports on A hy- 
peron production have been given elsewhere [12]. The S~ 
measurement reported here is unique in that it represents 
the lowest incident energy for which such a measurement 
has been made in heavy ion collisions (the laboratory 
energy threshold for direct production in NN collisions 
is 3.74 GeV). The production of multi-strange hadrons is 
arguably most interesting close to the production thresh- 
old due to the suppression of multistep processes (the 
S~ particle is mainly produced from the strangeness- 
exchange reactions {K{A,'E,) 2(77,77)). Collisions of 6 
AGeV Au + Au are expected to produce nuclear mat- 
ter with high baryon density and relatively low chemical 
freeze-out temperature. Thus, marked changes in the 
model predictions are expected for multi-strange particle 
yields. 

The experiment has been performed with the EOS 
TPC [13] at the Alternating Gradient Synchrotron (AGS) 
at Brookhaven National Laboratory (BNL). Details on 
the detector and its setup have been reported earlier 
[12-14] . The data presented here benefit from the excel- 
lent coverage, continuous 3D-tracking, and particle iden- 
tification capabilities of the TPC. These features are cru- 
cial for the efficient detection and reconstruction of A's 
and S~. 

Several studies have demonstrated the viability of ver- 
tex and trajectory reconstruction for V'^'s, i.e. from a 
track pair resulting from two charged daughter particles 
[12,15,16]. However, a signature for the doubly-strange 
S~ must be demonstrated by finding the two consecutive 
decays: S~ ^ Att^ followed by A ^ p7r~. This requires 
locating three correlated tracks and is considerably more 
difficult, especially for low yields near the threshold. The 
upper and lower panels of Fig. 1 show signature peaks in 
the invariant mass spectra for primary A's and for the 
S~ respectively. The mass peak is well centered at 
its rest mass of 1.32 GeV. The width of the distribution 
for E- (FWHM) is 10 MeV compared to 6 MeV for the 
width of the A as can be expected from the resolution of 
the TPC. 

A description for the reconstruction of primary A's is 
given in Refs. [12,15]. S~ reconstruction proceeded via 
a combination of two distinct steps. In the first, primary 
and secondary A-hyperons were reconstructed from the 
daughters of their charged particle decay, A — > p + n~ 
(branching ratio ~ 64%) following essentially the pro- 
cedure outlined in Refs. [12] and [15], except that the 
training (see below) was different for A's that do not orig- 
inate directly from the primary vertex (secondary A's). 
All TPC tracks in an event were reconstructed followed 
by the calculation of an overall event vertex. Thereafter, 
each p-jT~ pair was considered and their point of closest 
approach obtained. Pairs whose member trajectories in- 



tersect (with fairly loose criteria such as decay distance 
from the event vertex > 0.5 cm ) at a point other than 
the main event vertex were assigned as A candidates (pri- 
mary and secondary) and evaluated to yield an invariant 
mass and associated momentum. These A candidates 
were then passed to a fully connected feed-forward mul- 
tilayered neural network [16] trained to separate "true" 
secondary A's from the combinatoric background. The 
network was trained from a set consisting of "true" sec- 
ondary A's and a set consisting of a combinatoric back- 
ground. "True" secondary A's were generated by tagging 
and embedding simulated A's (created from the decay of 
S~ particles from RQMD [17] calculations) in raw data 
events in a detailed GEANT simulation of the TPC. The 
combinatoric background or "fake" A's were generated 
via a mixed event procedure in which the daughter par- 
ticles of the A {pn~) were chosen from different data 
events. 

In step 2, the S~ particle was reconstructed from its 
weak decay: E!~ — > Att" (branching ratio = 100%; cr = 
4.9 cm), using the secondary A's found in step 1. To do 
this, each Att" pair was evaluated to assign S~ candi- 
dates [and their associated invariant mass and momen- 
tum] which were passed to a fully connected feed for- 
ward multilayered neural network [16] trained to sepa- 
rate "true" S~'s from the combinatoric background. For 
this latter step, the network was trained from a set con- 
sisting of "true" S^'s and a set consisting of a combi- 
natoric background following a procedure similar to that 
described earlier for A's. 

For the four centrality selections (central to periph- 
eral) presented here 71656, 140706, 154542 and 57236 
events were processed to yield 78, 118, 100 and 16 S~'s 
respectively. The corresponding number of primary A's 
are 21526, 37090, 30973 and 6994 respectively. Rela- 
tively good phase space coverage was achieved for both 
the S^'s (0 < pT < 1-2 GeV/c, -0.5 < y„ < 0.5) and the 
A's (0 < pT < 1-2 GeV/c, -0.7 < ?y„ < 0.6). Here, y„ is 
the normalized rapidity. The detection efficiency for S~ 
ranged from 5.4 ± 0.4 x 10"''' to 23.5 ± 0.5 x 10"'''. Similar 
efficiencies for primary A's ranged from 2.3 ± 0.2 x 10"^ 
to 6.8 ±0.2 X 10-2. 

The neural network optimizes (during training) and 
then applies (during reconstruction) a multidimensional 
cut on topological variables or input neurons [16], to dis- 
tinguish true S-'s. In addition, several "hard cuts" (e.g. 
dcoA^TT, where dca is the distance of closest approach of 
a A, TT to the event vertex) were employed to investigate 
and minimize systematic uncertainties. Figure 1 shows 
the results obtained with such restrictions. In general, 
tighter "hard cuts" led to improvements in the signal-to- 
background ratio, but with significant reductions in the 
raw yield due to poorer detection efficiencies. Neverthe- 
less, the efficiency corrected yields showed essentially no 
dependence on these cuts. 

For the different cut conditions there are, of course, dif- 
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ferent efficiencies for track finding and thus for S~ and A 
reconstruction. These efficiencies have been determined 
by simulations that involve embedding E~ particles from 
RQMD calculations into events from the data sample and 
running this modified data set through the neural net- 
work. Subsequently the ratio of reconstructed S~ parti- 
cles in the real data divided by the simulated efficiency 
has been studied for a large variety of cut conditions. 
This method allows us to minimize systematic errors for 
the applied hard cut and thus we expect that statistical 
errors predominate. As a further test of this efficiency 
determination method, a model-independent integration 
over pt and rapidity was made for A's. The result was 
in very good agreement with the method used here for 
both A's and S~'s. 

Several centrality cuts have been made in order to ex- 
amine the variation of the S~ multiplicity on impact pa- 
rameter and the number of participating nucleons Apart- 
For each centrality cut, an impact parameter range was 
assigned via its respective fraction of the minimum bias 
cross section as described in Ref. [18]. The corresponding 
values for Apart were obtained via the Glauber model [19] . 
Figures 2a and b show the respective centrality (Apart) 
dependence of the multiplicity for S~ and (A + S*^) as 
the A yield includes feed-down contributions from the TP. 
The ratio S~/(A -|- is shown as a function of Apart 
in Fig. 2c. Figures 2a and b indicate a relatively strong 
increase in the S~ and (A -I- S") multiplicities with in- 
creasing Apart- However, this dependence is greater for 
the S~ as can be seen in Fig. 2c. An extrapolation [with 
the fit function multiplicity = const x Apart"] to the 
data shown in Fig. 2c leads to the ratio S~/(A -|- E°) 
of 0.017 ± 0.005 for most central collisions (6=0 fm or 

Apart ~ 370). 

Figure 2 also shows the results of calculations from 
the transport model RQMD [17] as indicated. For the 
single-strange hypcrons A + E*^, the model prediction 
of an increase in multiplicity with increasing centrality 
follows the trend of the data rather well. A similarly 
good agreement is apparent for the doubly-strange 
hyperon. This rather striking agreement between data 
and theory is in direct contrast to the observation that 
transport models are unable to describe the measured 
yields for S~ at SPS energies [11]. Next we compare the 
ratio S~/(A -(- E°) to that obtained at higher incident 
energies and to statistical model predictions as shown in 
Fig. 3. 

The statistical model makes predictions for this ratio 
as a function of ^/s [20]. The general freeze-out curve [21] 
taken together with the values of the chemical freeze-out 
temperature T and baryochemical potential /is as a func- 
tion of -^i [20] define all of the particle ratios. The dotted 
line in Fig. 3 indicates the ratio S~/(A -I- E") calculated 
for central collisions. It exhibits a pronounced rise in the 
AGS/SPS energy regime and then flattens off toward the 
RHIC domain of \/s > 100 GeV. Parameters for the sta- 



tistical model are determined from the systematics of fits 
to experimental data, and hence have finite uncertain- 
ties. Therefore, the curve in Fig. 3 should in principle 
be replaced by a band. Evaluation of these uncertainties 
as a function of energy is beyond the scope of this work. 

The rise up to SPS energies results from two primary 
effects. The first is an increase in temperature, from « 
50 MeV (SIS, GSl) up to « 170 MeV (SPS/RHIC). The 
second effect originates from the fact that one needs to 
use a canonical description at the lower incident ener- 
gies. For central collisions at SPS energies, the number 
of multi-strange particles is high and a grand-canonical 
description (global strangeness conservation) is sufficient. 
At lower incident energies, this number is much smaller 
requiring a canonical treatment (local strangeness con- 
servation) which leads to an additional suppression to- 
wards the lower incident energies [9]. Above top SPS 
energies the temperature is roughly constant, while 
decreases. Since the ratio S~/(A -I- E'^) is independent 
of j2B, it saturates (neglecting higher resonances which 
are taken into account in the model calculations) to the 
limit 5h-/(5a+.9s") ' exp[-{E^- + ns ~ E/^^j2o)/T] with 
gi the degeneracy, /xs the strangeness chemical potential 
and Ei the energy of particle i. Detailed calculations for 
the centrality dependence are not yet available. 

Data for E" production are indeed scarce. Results 
obtained at SPS from NA49 [22] and very recently at 
RHIC [23] for central coUisions of Pb+Pb/Au+Au are 
shown as filled symbols in Fig. 3. The RHIC point [for 
midrapidity only] agrees (well) with the statistical model. 
The SPS point for the 4tt integrated yield is below the 
statistical-model curve. At the highest AGS energy a 
data point (open circle) is shown for a Si-|-Au measure- 
ment performed at slightly forward rapidities (y =1.4- 
2.9) [24]. This data point exceeds the curve, but is within 
two standard deviations even neglecting the systematic 
uncertainty in the statistical model calculations. The 
filled star in Fig. 3 indicates our new Air integrated result 
for central Au-|-Au collisions at 6 AGeV. This value is 
in good agreement with the statistical model prediction 

[20]- .... 

This data point is also in good agreement with the 
predictions of the RQMD model (cf. Fig. 2c). Thus, 
we conjecture that hadronic processes are able to explain 
the measured yields. This is in agreement with recent 
UrQMD calculations [25] which give a detailed account 
of the channels responsible for 5~ production. Since the 
statistical model predictions also give a good represen- 
tation of the data, the question is raised as to how far 
transport model calculations are from chemical equilib- 
rium. Additional theoretical work is clearly needed to 
answer this question and to unravel the production mech- 
anism/s for multistrange hadrons from energies close to 
threshold up to the SPS [where current calculations fail 
to reproduce the measured yields by hadronic processes 
alone] . 
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In conclusion, we have presented the first measurement 
for the doubly strange hyperon S~ in near-threshold 
Au+Au collisions. The multiplicity of the S~ is observed 
to increase more rapidly with centrality than the multi- 
plicity of A -|- but the measured ratio S^/(A -|- E°) 
for central collisions agrees well with both statistical and 
transport model predictions. This agreement may be an 
indication that despite the rather short reaction times at 
AGS energies, the available phase space for S~ is essen- 
tially filled. 

This work was supported in part by the U.S. De- 
partment of Energy under grants DE-FG02-87ER40331, 
DE-FG02-89ER40531, DE-FG02-88ER40408, DE-FG02- 
87ER40324, and contract DE-AC03-76SF00098; by the 
US National Science Foundation under Grants PHY-98- 
04672, PHY-9722653, PHY-96-05207,PHY-9601271, and 
INT-9225096; by the University of Auckland Research 
Committee, NZ/USA Cooperative Science Programme 
CSP 95/33; and by the National Natural Science Foun- 
dation of P.R. China under grant 19875012. 



(1997) . 

[17] H. Sorgo, Phys. Rev. C52, 3291 (1995). 

[18] P. Chung ct al., Phys. Rev. C66, 021901(R) (2002). 

[19] R. J. Glauber, Phys. Rev. 100, 242 (1955); 

[20] P. Braun-Munzinger, J. Cleymans, H. Oeschler, K. Red- 

lich, Nucl. Phys. A697, 902 (2002). 
[21] J. Cleymans and K. Rodlich, Phys. Rev. Lett. 81, 5284, 

(1998) ; Phys. Rev. C60, 054908 (1999). 

[22] Latest results with 4n yields are from: S.V. Afanasiev et 

al., (NA49 Collaboration), Phys. Lett. B538, 275 (2002); 

A. Mischke ot al., (NA49 Collaboration), (SQM2001), J. 

Phys. G: Nucl. Phys. 28, 1761 (2002). 
[23] J. Castillo et al., (STAR Collaboration), Quark Matter 

2002. 

[24] S.E. Eiseman et al., (E810 CoUab.), Phys. Lett. B325, 
322, (1994). 

[25] M. Reiter et al., talk presented at the Trento workshop on 
transport model calculations, Trento, May 19th - 23rd, 
2003 and to be published. 



[1] J. Rafelski, Phys. Lett. B97, 297 (1980); 

P. Koch B. Miiller and J. Rafelski, Phys. Rep. 142, 168 
(1986). 

[2] J. Rafelski, Phys. Lett. B262, 333 (1991). 

[3] Proc. of the 6th International conference on Strange 

Quarks in Matter, J. Phys. G: Nucl. Part. Phys. 28, 1517 

(2002). 

[4] E. Andersen et al., (WA97 Collaboration), Phys. Lett. 
B449, 401 (1999) ; N. Carrer, (NA57 Collaboration), 
QM2001. 

[5] P. Braun-Munzinger, J. Stachel, J. P. Wessels and N. 
Xu, Phys. Lett. B344, 43 (1995) and B365, 1 (1996); P. 
Braun-Munzinger and J. Stachel, Nucl. Phys. A606, 320 

(1996). 

[6] P. Braun-Munzinger, I. Heppe, J. Stachel, Phys. Lett. 

B344, 43 (1995). 
[7] P. Braun-Munzinger, D.J. Magestro, K. Redlich and J. 

Stachel, Phys. Lett. B518, 41 (2001). 
[8] e.g. International conference on "Quark Matter 2002", 

Nantes, France, July 2002. 
[9] J. Cleymans, H. Oeschler and K. Redlich, Phys. Rev. 
C59, 1663 (1999); Phys. Lett. B485, 27 (2001). 
[10] J. S. Hamieh, K. Redhch and A. Tounsi, Phys. 

Lett. B486, 61 (2000). 
[11] S. Soff ct al., Phys. Lett. B471, 89 (1999). 
[12] P. Chung et al., J. Phys. G 25, 255 (1999); C. Pinkenburg 
et al, Proc. of Quark Matter 2001, Stony Brook, New 
York, January 2001, Nucl. Phys. A 698, 495c (2002). 
[13] G. Rai ot al., IEEE Trans. Nucl. Sci. 37, 56 (1990). 
[14] C. Pinkenburg et al., Phys. Rev. Lett. 83, 1295 (1999). 

(1999); P. Chung et al., Phys. Rev. Lett. 85, 940 (2000). 
[15] M. Justice, et al., Phys. Lett. B440, 12, (1998). 
[16] M. Justice, Nucl. Instrum. Meth. Phys. Res. A400, 463 



Au(6AGeV) + Au(b<7fm) 



(a) A 




1.29 1.3 1.31 1.32 1.33 1.34 1.35 



Inv. Mass (GeV/c^) 



4 



FIG. 1. Invariant mass spectrum for (a) A and (b) S~ in 6 
AGeV semi-central Au+Au collisions {b < 7 fm). The curve 
in panel (b) is a Gaussian fit to the data points. The Full 
Widths at Half Maximum for the peaks are indicated. With 
these cuts there are ~ 250 E~ counts (net of background). 



FIG. 2. Panels (a) and (b) show the Apart dependence of 
the multiplicity of H~ (circles) and (A + E'') (squares), respec- 
tively, from 6 AGeV An -I- Au collisions. Central collision have 
Apart ~ 370. Panel (c) shows the corresponding multiplicity 
ratio (stars), H~/(A + S°), as a function of Apart. Values from 
RQMD calculations are indicated via the dotted lines. 



0.3 
0.2 
0.1 

1 

1 10 
5 


g 0.02 

CO 

^ 0.01 



_afiii Sfm 4firi ZfuL 



RQMD (V2.3) 



+ 



(b) 



(A + E° 



(c) 



★ E7(A + r 



50 100 150 200 250 300 350 400 
^part 



10 



+ 



10 



< 


, , ^ 

S RHIC 




SPS 


AGS 


I 


★ E895 Au(6 AGeV)+Au 
O E810Si(14.6AGeV)+Au 

NA49 Pb(158 AGeV)+Pb 
▲ STAR Au+Au 

Statistical Model 



10 

Vs (GeV) 



10 



FIG. 3. The multiplicity ratio of H"/(A -I- E°) as a func- 
tion of -^/(s) Filled symbols correspond to central collisions of 
Au+Au/Pb-|-Pb. The datum indicated by the open circle is 
for the asymmetric system Si-|-Au. The dotted line shows re- 
sults from the statistical model for central Au-I-Au collisions. 
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